Constraining Lorentz violation with cosmology 
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The Einstein-Aether theory provides a simple, dynamical mechanism for breaking Lorentz in- 
variance. It does so within a generally covariant context and may emerge from quantum effects in 
more fundamental theories. The theory leads to a preferred frame and can have distinct experimen- 
tal signatures. In this letter, we perforin a comprehensive study of the cosmological effects of the 
Einstein- Aether theory and use observational data to constrain it. Allied to previously determined 
consistency and experimental constraints, we find that an Einstein- Aether universe can fit experi- 
mental data over a wide range of its parameter space, but requires a specific rescaling of the other 
cosmological densities. 



I. INTRODUCTION 

The spacetime symmetry of local Lorentz invariance is 
a cornerstone of modern physics [l|, but is not inviolate. 
Violations can occur in quantum gravity theories, with 
the symmetry emergent and approximate at macroscopic 
levels 0- In the particle physics sector the symmetry has 
been experimentally verified to extremely high precision 
Q . On the large scales characteristic of the gravitational 
sector, however, constraints are much less certain. In this 
letter we explore the extent to which precision cosmology 
can constrain a Lorentz-violating theory. 

The theoretical workhorse for studying violation of 
Lorentz symmetry in gravitation is the Einstein- Aether 
theory [4j, a simple, elegant proposal for dynamically 
violating Lorentz invariance within the framework of a 
diffeomorphism-invariant theory. It is a refinement of 
the gravitationally coupled vector field theories first pro- 
posed by Will and Nordvedt in 1972 [5] and has been 
explored in exquisite detail by Jacobson, Mattingly, Fos- 
ter and collaborators fg, 7, 8]. A Lorentz-violating vector 
field, henceforth called the aether, will affect cosmology: 
it can lead to a renormalization of the Newton constant 
Q, leave an imprint on perturbations in the early uni- 
verse [l3, El , and in more elaborate actions it may even 
affect the growth rate of stucture [H, [l3| . Preliminary 
calculations have been done on the affect of the aether on 
anisotropics of the cosmic microwave background [14j . 

The action for the Einstein- Aether is: 
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where g ab is the metric, R is the Ricci scalar of that met- 
ric, Sm is the matter action, and £ is constructed to be 
generally covariant and local. G is the bare gravitational 
constant, not necessarily equal to the locally measured 
value. Sm couples only to the metric g a b and not to A b 
and 
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|15j |. We will use the notation c%2... = c\ + c-i + ... The 
gravitational field equations for this model take the form: 
G a b = T ao + 87rGT a (, where the stress-energy tensor for 
the vector field T a b is given in [4[ and T ab describes the 
conventional fluids. 

A number of constraints on the CjS have been de- 
rived. Most notably a Parametrized Post Newtonian 
(PPN) analysis of the theory leads to a reduction in the 
dimensionality of parameter space such that C2 and C4 
can be expressed in terms of the other two parameters: 
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ditionally, the squared speeds of the gravitational and 
aether waves with respect to the preferred frame must 
be greater than one so as to prevent the generation of 
vacuum Cerenkov radiation by cosmic rays [161 ] . We shall 
label this space of models as C. A final constraint arises 
from considering the effects of the aether on the damp- 
ing rate of binary pulsars. The rate of energy loss in 
such systems by gravitational radiation agrees with the 
prediction of General Relativity to one part in 10 3 . It 
has been shown [l7| that, for the Einstein- Aether theory 
to agree with General Relativity for these systems, we 
require that c + = c\ + C3 and c_ = c\ — C3 arc related 
by an algebraic constraint (shown as the dashed line in 
figure [3]) [HI]. A more exotic, but viable, subset of the 
parameter space can be considered in which c\ = C3 = 0. 
The PPN and pulsar constraints do not apply here and a 
cosmological analysis is potentially the only way of con- 
straining the values of the coupling constants. We shall 
label this alternate space of models as £. Note that, 
in what follows, we will write down the equations in a 
general form and then study the two subspaces C and £ 
independently. 

We now focus on cosmological scales and assume 
a homogeneous and isotropic background spacetime in 
which the metric is of the form g ab dx a dx b = —dt 2 + 
a{t) 2r )ijdx l dx : ' where t is physical time, 7^ is the iden- 
tity matrix, and a(t) is the scale factor. Throughout this 
letter, subscripts i and j will run 1 to 3. The vector 
field must respect the spatial homogeneity and isotropy 
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of the system and so will only have a non-vanishing T 
component; this constraint fixes A b = (1,0,0,0). The 
energy-momentum tensor of the matter will include the 
standard menagerie: photons, neutrinos, baryons, dark 
matter, and the cosmological constant. 

In this background, the t-t component of the aether 
stress energy tensor is equal to (3/2)aH 2 where a = c\ + 
3c2 + C3. For models C we have a — — 2c + c_/(c + + c_). 
The fractional energy densities in the various compo- 
nents are given by f2j(a) = 8itG pi(a) / ZHq and VI ae = 
(a/2)Eif2i/(l - a/2) where H = 100ft km s^Mpc" 1 
is the Hubble constant today and pi is the energy den- 
sity in the fluid component i. In quasistatic space- 
times the aether exhibits tracking behaviour such that 
the locally measured value of Newton's constant is ac- 
tually G N = G/(l + c-u/2) 0. For models C we have 
C14 = 2c+c_/(c+ + c_) [20|. Thus, given a value pi, the 
actual is related to the value Clm inferred using a 
locally measured value of Gn as: Qj = (1 + Ci4/2)f2jvi- 
Hence, explicitly using our expression for fi^g, the Fried- 
mann equation becomes: 
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To fully explore the cosmological consequences of the 
aether, we must consider linear perturbations around the 
background. We will do so in the synchronous gauge 
and use conformal time coordinates: given g a bdx a dx b = 
—a 2 dr 2 + a 2 [7y + h i j]dx l dx : > , the two scalar potentials 77 

and h are defined by: ft.y(x, r) = d 3 fce x [fcjfcj/i(k, r) + 
(kikj — i<5y)6?7(k, r)]. The aether field can be written as 
A d = -(1, 9jV); the zeroth component is fixed equal to 
a -1 by the gauge choice and the fixed-norm constraint. 
Instead of V itself, we choose to to use the variable: £ = 
V — (h + 6r))' with which the evolution equations take 
a particularly instructive form. 

The gravitational field equations are 

(1 - ^a)k 2 i = 4irGa 2 iki6T° 3 + ifc 4 c 123 £ 

and 
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For models C we have C123 = 20^/3(0+ + c_). 
The aether equation of motion is: 
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where H is the conformal Hubble parameter, primes are 
derivatives with respect to r and E/ = —8irGa 2 (kik^ — 
|^" 7 i )E t ^, where E^ is the traceless component of the fluid 
stress energy tensor. The homogeneous 'sourceless' solu- 
tion to the above equation during an era where a cx r™ 
is t(k,r) = T^iM^J^Cskr) + f 2 {k)Y(p,c s kT)] 
where are functions to be fixed by boundary condi- 
tions, J and Y are Bessel functions and the various con- 
stants are defined through: b\ = —2n — (cu + a)(n 2 + 
n)/[c 14 (l + c+)], = (1 - 8n + b 1 + 6 2 ) 1 / 2 an d c 2 = 
c i23/(ci4(l + c + )). With c 2 positive, the solutions are 
damped and oscillatory solutions when c s kr 3> 1 and 
power law when c s kr <SC 1 [l2l ]. 

It was shown in [ici| that the primordial scalar power 
spectrum T-^ (where $ is the trace perturbation to the 
metric in the conformal Newtonian gauge) is modified rel- 
ative to that of a Universe with no aether, P$, through 
"i-— c+1 2 

1 a . 14 - — whilst £ and £' are driven to a van- 

l+c+ J s ' 

ishingly small value compared to their values at the onset 
of inflation. We will work with the equivalent initial con- 
ditions in the synchronous gauge. 

To study these effects in detail, we have modified the 
Boltzmann code CMBEASY (2lj by adding a Newton- 
Raphson solver for the Hubble parameter, and including 
the aether components in the density and pressure; the 
perturbation evolution has been modified by adding £ 
and £' as the integrated components. In Figure[T]we show 
the effect of the aether on the angular power spectrum 
of anisotropies in the CMB and the power spectrum of 
galaxies (we superpose the WMAP and SDSS data) for 
a selection of parameters in class C. 

The dominant effect for smaller values of c + is on the 
large-scale CMB, through the integrated Sachs- Wolfe ef- 
fect; it leads to a supression on large scales (which curi- 
ously enough is favoured by large scale CMB data). As 
expected, the overall growing mode of matter perturba- 
tions is very weakly affected and the change on the power 
spectrum of galaxies is marginal. 

As is usual in cosmological model testing, we com- 
pute parameter constraints using a Monte-Carlo Markov 
Chain (MCMC)[23j. We explore a 6 dimensional pa- 
rameter space consisting of the fractional baryon den- 
sity, Of,, the fractional matter density, fijvf, the Hubble 
constant, Hq, the scalar spectral index, n$, the optical 
depth, td, the overall amplitude of fluctuations, the bias 
factor of SDSS galaxies and the two aether parameters, 
c + and c_. We constrain parameters using the WMAP 
3-year release, the Boomerang 03 release and data from 
ACBAR and VSA[H, as well as the SDSS and 2DF 
surveys [IjJ [26| . We also use measurements of the lumi- 
nosity distance as a function of rcdshift from supernovae 
la measurements [27| but have found that these data sets 
have very little ability to constrain this class of models. 
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FIG. 1: The angular power spectrum of the CMB (bottom) 
and the power spectrum of galaxies (top) for a sample of class 
C Einstein- Aether models, with different c+ (with c_ chosen 
to satisfy the weak field binary pulsar constraint — the dashed 
line of figure [4} The other parameters have their AC DM best 
fit values, with the fii rescaled as described in the text. Su- 
perposed are the WMAP and SDSS datasets. 
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FIG. 2: The angular power spectrum of the CMB (bottom) 
and the power spectrum of galaxies (top) for a sample of ex- 
otic class £ Einstein- Aether models where ci = C2 = C3 = 0. 
The other parameters have their ACDM best fit values, with 
the radiation density modified to account for the change in 
the gravitational constant. Superposed are the WMAP and 
SDSS datasets. 

The marginalized constraints from the CMB and large- 
scale structure on the two aethereal parameters in model 
C are shown in hgure [4] [lj| . The best-fit aether model is 
mildly superior to standard ACDM cosmology, at about 
2<t, at a cost of two extra parameters. 

The soft lower limit at c_ > —0.5 comes from a prior 
on the baryon fraction. This signals an important char- 
acteristic of these models: the strong correlation between 
the fractional energy density in the aether, Oab, and the 
other energy components. This is perhaps the primary 
result of our analysis and is illustrated in figure [3J the 
CMB and LSS data restrict the background to evolve as 
in the ACDM case, which in turn leads to a rescaling of 
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FIG. 3: Joint constraints on the fractional aether density, 
Qae , with the physical dark matter density, 0, c h 2 , the phys- 
ical baryon density, fij/i 2 and the fractional A density, Oa- 
The contours are 1 and 2a. 

the different energy components in the presence of the 
aether. Naturally this also affects the constraints on the 
other cosmological parameters. These constraints, un- 
der the ACDM and aether models with and without the 
weak binary pulsar constraint, are shown in table [T] As 
expected, the largest shift is seen in the Hi. 
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FIG. 4: Joint constraints on the parameters — c+ and — c_. 
The black lines are the 1 and 2a contours, where we have 
marginalized over the values of the parameters. The hatched 
region is excluded by Cerenkov constraints; the dashed line 
indicates where weak-field constraints from binary pulsars are 
met. Both are taken from Q. 

As stated above, the CMB and LSS play the dominant 
role in generating these constraints, and interestingly 
enough this is through the change in the background evo- 
lution and its effect on the metric perturbations, and not 
necessarily through the presence of perturbations in the 
vector field. Indeed, artificially switcing off the pertur- 
bations in the aether field has essentially no effect on 
the power spectrum of LSS and a small effect (of ap- 
proximately 10%) on the angular power spectrum of the 
CMB. 

So far we have focused on models in class C, where 
we found that the coupling constants are allowed to vary 
quite widely. In the case of models in class £ , the cosmo- 
logical data are far more restrictive. For example, fixing 
Ci = C2 = C3 = we find that —0.05 < C4 < (note 
that in this case £Iae = 0)- If we allow ci to be non-zero 
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Parameter 


ACDM 


General 


Weak Pulsar 


n c h 2 


0.137 ±0.004 


0.097 ±0.01 


0.098 ±0.01 


Qbh 2 


0.022 ± 0.001 


0.019 ±0.002 


0.019 ±0.002 


H 


69.7 ± 1.6 


71.7 ±2.0 


72.5 ±2.4 


TD 


0.08 ±0.029 


0.077 ±0.027 


0.078 ± 0.028 


n s 


0.955 ±0.015 


0.976 ±0.02 


0.984 ± 0.024 


Qa 


0.671 ±0.019 


0.61 ± 0.05 


0.67 ±0.028 


Cl 




-0.46 ±0.14 


-0.26 ± 0.12 


C2 




0.34 ±0.1 


0.20 ± 0.09 


C3 




-0.23 ±0.1 


-0.12 ±0.05 


C4 




0.13 ±0.09 


0.05 ± 0.02 



TABLE I: Mean and la error values of the marginalized like- 
lihoods for a range of cosmological parameters. The left most 
column is for ACDM with no aether, the central column for 
general class C models, and the right hand column for class 
C models with the weak pulsar constraint (on the dashed line 
of Figure 1). 

we find that both c-i and C4 must be in [—0.01,0]. The 
reason for this constraint is illustrated in figured) at low 
£ the integrated Sachs- Wolfe effect induced by the mod- 
ified potentials is large enough to disrupt the Ci. These 
are the strongest constraints on these parameters that 
currently exist. 

In this letter we have studied the effect of Lorentz vi- 
olation on cosmology as parametrized by the Einstein- 
Aether model. We have found that Lorentz violation in 



this form is compatible with current cosmological data 
and, combined with other non-cosmological probes we 
have found constraints on c+ and c_. The data also 
require the rescaled combination of density parameters, 
in which the background evolution is unchanged from a 
ACDM universe. We have also found tight constraints 
on the other allowed range of parameter space, £, which 
has, until now, been relatively unconstrained by other 
methods. Collectively these constraints arise from tests 
on distance scales spanning more than fifteen orders of 
magnitude. 

There are of, course, other possible ways of param- 
eterizing Lorentz violation which are not encompassed 
by the Einstein-Aether model. In particular one may 
relax the fixed-norm constraint on the aether field [28| 
or couple it directly to the matter content of the Uni- 
verse [3|. Such theories tend to have a much stronger 
effect on the evolution of the background cosmology [29| 
or lead to distinct experimental signatures [30(. Hence 
the results presented here are currently the most com- 
prehensive (though conservative) constraints on the CjS, 
and thus on Lorentz- violating vector theories. 
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